Strong impact of light induced conical intersections 
on the spectrum of diatomic molecules 
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We show that dressing of diatomic molecules by running laser waves gives rise to conical inter- 
sections (CIs). Due to presence of such CIs, the rovibronic molecular motions are strongly coupled. 
A pronounced impact of the CI on the spectrum of Na2 molecule is demonstrated via numerical 
calculation for weak and moderate laser intensity, and an experiment is suggested on this basis. The 
position of the light induced CI and the strength of its non-adiabatic couplings can be chosen by 
changing the frequency and intensity of the used running laser wave. This offers new possibilities 
to control the photo-induced rovibronic molecular dynamics. 
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The spectrum of diatomics in weak laser fields is 
usually interpreted in terms of the Franck-Condon 
(FC) principle (see e.g. Ref. [lj]). If the laser light be- 
comes more intense, the common approach is to take 
into consideration the two relevant electronic poten- 
tial energy curves (PECs) coupled by the laser, and to 
account for rotational transitions via the selection rule 
AJ = ±1 (see e.g. Refs.JH-Ql and citations therein). 
For even stronger laser intensities, molecular rotations 
must be accounted for, as done e.g. in Ref. 

For standing laser waves we have shown [fl that 
laser light induces conical intersections (CIs) which 
couple the center of mass motion with the internal 
rovibronic degrees of freedom. As we shall discuss in 
the present work, laser induced CIs are formed even 
for running laser waves. It is well known Q that for 
the formation of a CI one needs at least two nuclear 
degrees of freedom whose changes affect the electronic 
wavefunction. These are not available for a free di- 
atomic molecule. However, as we shall demonstrate 
below, the rotation becomes the missing degree of free- 
dom to allow for the formation of CIs for diatomics in 
the presence of a laser field. Since Na,2 is a favorite 
and well studied molecule of laser physics @, we 
choose it as our example. The situation is best illus- 
trated in Fig.l. The top panel shows the potential en- 
ergy curves of the X 1 S+ and A *E+ electronic states 
of Na,2 ■ The vertical arrow labelled by hu> l indicates 
the transition between these two electronic states by 
a weak laser field in accordance with the FC theory. 
Fig. lb depicts the adiabatic potentials obtained by 
dressing the just mentioned two PECs by the laser 
following the standard approach of Rcfs.JJQ referred 
to above. In Fig.lc we plot the laser induced CI re- 
sulting from our approach as discussed below. 

The presence of CIs is well known in polyatomics to 
exert an enormous impact on the molecular dynamics 
@, E3] ■ We will show here for the first time that this 



is also the case for the CI induced by running laser 
waves in diatomics. 

Let us briefly describe now the underlying theoret- 
ical background. Using the same formalism and as- 
sumptions as in Q one can derive the following Hamil- 
tonian for rovibronic molecular motions: 
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Here, R and {6, if) are the molecular vibrational and 
rotational coordinates, respectively, fx stands for the 
reduced mass of the diatomic, and Lg^ represents the 
angular momentum operator of the nuclei. The cen- 
ter of mass motion has been eliminated via the usual 
adiabatic separation approach which is physically ade- 
quate for the case of running waves. Subscripts X and 
A refer to the two electronic states coupled by the laser 
(whose frequency is lot, and amplitude £q). To be spe- 
cific we assume that these two electronic states possess 
E type symmetry, such that the molecule-light cou- 
pling term is proportional to d(R) cos 9, where d(R) 
is the corresponding dipole transition matrix element 
in the body fixed frame. As discussed in detail in our 
previous work [(|, the potential matrix of the Hamil- 
tonian ([TJ gives rise to a conical intersection whenever 
the conditions Vx(R) = Va(R) — Tiul and 8 — ir/2 are 
met. This situation is shown explicitly in Fig.lc. 

The conventional Hamiltonian H widely used for 
diatomics in laser fields (see e.g. Refs.0-01) is obtained 
from ([!} by replacing cos 8 by its optical transition 
matrix element 
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Here, \l) = Pi(cos6) are the Legendre polynomials. 
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FIG. 1: Electronic potential energies of the Na,2 molecule 
dressed by laser light, pertaining to different levels of the- 
ory discussed in the text. Top: The field free PECs. Mid- 
dle: The two lowest lying dressed adiabatic potentials for 
a field intensity of 3 • 10 s W/cm 2 . Bottom: Conical inter- 
section induced by a running laser wave. For visualization 
reasons we have used here a rather strong laser intensity of 
10 14 W/cm 2 . Note however that the conical intersection 
exists for any nonzero field strength. 

To illustrate the impact of the light induced CI on 
observable quantities we discuss the Nd2 molecule, 
considered above in Fig.l. The PECs and the tran- 
sition dipole functions pertaining to the^X 1 E+ and 



121 and the 



A 1 E+ electronic states are taken from 
laser wavelength is A = 663 nm. 

Solution of the eigenvalue problem of the Hamil- 
tonian ^ provides the quasi-energy eigenvalue spec- 
trum. These quasi-energy eigenvalues and eigenfunc- 
tions are calculated numerically by representing the H 
as a matrix in the basis set of the free molecular rovi- 
bronic states. An energy truncation criterion is used 
to reduce the number of basis functions to a manage- 
able size. Numerical convergence with respect to this 
energy cutoff has been checked and achieved. We shall 
thus hereafter refer to the obtained numerical results 
as " exact" . 

In the case of weak laser intensities one expects that 
the second order perturbation theory (PT2) in £q will 



be adequate for the energy levels. Indeed, if we use 
the common approach (which we call in the following 
the "no CI" case), we find this expectation to be ful- 
filled. The situation changes dramatically when full 
account of the CI is taken. This is apparent in Fig. 2, 
where the deviations of the PT2 results from the ex- 
act ones are shown. Breakdown of PT2 in the vicinity 
and above the energetic location of the CI is clearly 
demonstrated. 
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FIG. 2: Breakdown of the second order perturbation the- 
ory in the vicinity of the energy of the CI and above it. 
The position of the CI is marked by a vertical dotted line. 

In Fig. 2 all the lowest energy levels of the Hamilto- 
nian (JT|) have been considered. To make contact with 
a possible experiment, we study below the absorption 
spectrum of Na,2 where the levels are weighted by 
their intensities. We discuss three types of spectra: 
The FC spectrum obtained by employing the stan- 
dard FC principle, the "no CI" spectrum determined 
by the standard approach, and the "exact" spectrum 
following from ([1} . The Golden rule expression for the 
"exact" spectral intensity is given by \p(E)\ 2 , where 

p(E) = (*$(R,6) \d(R) cos 6 \*i(R, 6)) 

+ (*$(R,0)\d(R) cos (R, 0)) . (3) 

Here, $o(R,9) = (V$(R,6), *o (#> d )) T stands for 
that particular eigenstate of the Hamiltonian (Q} 
which reduces to the ground state of the free diatomic 
molecule in the limit of £q — > 0. Similarly, ^e(R, 9) = 
(^^(R, 6),^g{R, 9)) T represents any other eigenstate 
of the Hamiltonian ([1}, with quasi-energy E. The 
physical picture behind the just defined spectral in- 
tensity is as follows. First, a sample of free diatomic 
molecules (cold enough to be predominantly in their 
ground states) is dressed by an adiabatically switched 
on CW laser pulse with central frequency wj,. This 
leads via an adiabatic passage to the dressed molecular 
state ^o(R, 9) described above. Subsequently, the sys- 
tem is probed by another laser which facilitates dipole 
transitions from 1 i l Q(R, 9) to ^e(R, 9) according to the 
standard Golden rule theory as described by Eq. ([3]). 
The FC spectrum is obtained from ([3]) by identify- 
ing the initial state ^q(R,9) with the ground state 
of a free diatomic molecule, and by choosing the fi- 
nal states ^£;(i?, 9) to be the eigenstates of the A 1 E+ 
state in absence of the laser. It is known [l[ that the 
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FC spectrum becomes precise in the weak field limit 
of £q 0. For the "no CI" case, the spectrum is 
obtained from ([3]) by substituting for ^ (R,9) and 
$?e(R, 0) the corresponding eigenstates of the modi- 
fied approximative Hamiltonian H which was defined 
above using Eq. J2]). 
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FIG. 3: Rovibronic absorption spectra of Na,2 dressed 
by a running laser wave (A = 663 nm, intensity 3 • 10 8 
W/cm 2 ). The calculations are performed at different lev- 
els of theory (FC, "no CI", exact) discussed in the text. 
Top: overview of the spectrum. Bottom: Zoom into the 
range of 15400 and 16200 cm" 1 (with an inset showing 
the range between 15900 and 16200 cm - in an extended 
scale) . Fingerprints of CI are clearly demonstrated by pro- 
nounced deviation of the exact spectrum from the FC and 
"no CI" approximations. 

The computed spectra are depicted in Fig. 3 for a 
weak laser intensity of 3 • 10 8 W/cm 2 . Fig. 3a gives 
an overview of the calculated transition intensities. 
For most lines the FC approximation describes the 
spectrum qualitatively well. Nevertheless, clear devi- 
ations from the "no CI" spectrum are seen for most 
lines, with the intensity of the FC lines being higher 
(in particular for lines with energies about 18000 and 
28000 cm" 1 above the ground state). Moreover, many 
FC lines are seen to split into two in the "no CI" spec- 
trum. Interestingly, even at this weak laser field one 
can see obvious differences between the "no CI" and 
the " exact" spectra. The richness of the " exact" spec- 
trum becomes apparent on an expanded energy scale. 
In Fig. 3b we show the vicinity of the FC line between 
15400 and 16200 cm" 1 . While there is only a single 
FC line (characterized by the electronic state A, and 
u = 6, J = 1), we observe two lines in the "no CI" 



spectrum. The second weaker line is mainly the vi- 
brational level (y = 4, J = 0) of the lower electronic 
state X and has acquired its intensity by borrowing 
it from the more intense line {y = 6, J = 1) of the A 
electronic state. This borrowing of intensity is a well 
known non-adiabatic effect [Io|. Interestingly, the CI 
gives rise to further substantial splittings. Each of the 
"no CI" peaks splits into essentially three lines. For 
the weak laser intensity at hand, the splitting is mostly 
due to rotations. The assignment of the exact levels 
in Fig. 3b refers to the leading configuration in the ex- 
pansion of the corresponding eigenstates in terms of 
the free molecular states. In reality, several (or even 
many) configurations contribute. For example, the 
state associated with the most intense borrowing line 
in the spectrum (labelled by X, v = 4, J = 4) has a 
weight on the excited A electronic state of 20%. The 
averaged rotational quantum number Jx on the X 
electronic state equals to Jx ± A Jx = 1.8 ± 1.8 where 
AJx is the standard deviation. Analogously, on the 
A electronic state, J a ± AJ^ = 0.4 ± 1.2. 

We note that the absorption profiles depicted in 
Fig. 3 require spectroscopic resolution better than 50 
cm" 1 which can be obtained in the lab (see, e.g., 
Ref. @). In supersonic beam experiments Q the 
predicted phenomenon should be possible to observe. 
Fig. 3 may thus serve as a proposal of a future ex- 
periment aimed at demonstrating the impact of light 
induced CIs on molecular rovibronic spectra even for 
rather weak fields. 

Let us discuss now the situation encountered at the 
larger field intensity of 10 11 W/cm 2 . Here, all the 
molecular degrees of freedom are strongly coupled, 
and the corresponding eigenstates of the Hamiltonian 
(p} cannot be assigned to their field free counterparts. 
This complicates the use of the Golden rule formula 
([3]), as we are unable to populate the initial state 
^>o(R,9) via an adiabatic passage. In order to gain 
an insight into the nature of the obtained dressed 
molecular states, we instead consider the projection 
P(E) = |(*f (R,9)\(t>x =0 j =0 {R))\ 2 of each dressed 
state ^e(R, 0) on the field free molecular ground state 
0^ =OJ=o (i?). The quantity P(E) possesses a clear 
physical meaning. It is the probability of populat- 
ing the dressed state ^e{R, 0) when the laser pulse 
is imposed abruptly on the sample of free diatomic 
molecules residing in their ground states 4>v=o.J=o(R)- 
For the considered field intensity, pulse shaping tech- 
niques allow the production of nearly any square laser 
pulse, depending on the availability of a source that 
can cover the required bandwidth [ll| . 

Fig. 4 shows the obtained profile of P{E). Surpris- 
ingly, it turns out that the " no CI" spectrum has no re- 
semblance whatsoever with the exact spectrum. The 
" no CI" spectrum has become very sparse and the in- 
tensity of its lines is typically more than one order 
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FIG. 4: Probability P(E) of populating different dressed 
molecular states when an intense laser light (A = 663 nm, 
intensity 10 11 W/cm 2 ) arrives abruptly to a sample of free 
Na,2 molecules in their ground states. For visualization 
reasons, the exact spectrum is multiplied by factor 5. An 
inset shows a zoom into the region between 20000 and 
40000 cm -1 , where a completely irregular pattern of the 
profile of P(E) is evident. 

of magnitude larger than that of the lines of the ex- 
tremely dense exact spectrum. This sparseness is a 
strong indication that at this laser field strength the 
non-adiabatic effects in the "no CI" model have re- 
duced enormously. The high density of lines in the 
exact spectrum, on the other hand, is a signature 
of the sin gula r non-adiabatic coupling induced by a 
strong CI [10j. In particular, the spectrum above the 
energy of about 20000 cm -1 exhibits a typical irreg- 
ular pattern like the spectra of CIs in polyatomics 
which have been shown to behave like the spectra 
of random matrix ensembles [l3|.This irregularity re- 
flects the complete mixing of multiple configurations 
in the respective eigenstates. Our numerical analy- 
sis shows that both the X and A electronic states 
are essentially equally populated in all the eigenstates 
above the CI energy. Similarly, also the field free vi- 
brational as well as rotational levels are enormously 
mixed. For example, for the lines between 40000 and 
45000 cm" 1 one typically finds J x ± AJ X = 20 ± 20 
and J A ± AJ A = 25 ± 25, while v x ± Ais x = 3 ± 3 
and vx ± Ai/x — 5 ± 5. In all cases, the magnitudes 
of the calculated averages are roughly equal to their 
standard deviations. 

We have shown that even diatomic molecules ex- 
hibit CIs. These are induced by running laser waves, 
and live in the space of the molecular rotational and 
vibrational degrees of freedom. We have demon- 
strated that these CIs have a strong impact on the 
molecular spectrum even for weak laser fields. For in- 
termediate laser fields this impact becomes dramatic 
and the resulting spectrum has no resemblance with 
those computed by standard approaches. The cou- 
pling of two electronic states via the light induced 
conical intersection is substantially and sometimes 
even dramatically stronger than without it, and con- 
sequently, one can apply rather weak fields to achieve 
strong coupling between these two states. If needed, 



e.g., when the field is strong and other electronic states 
are close by, one can easily include more electronic 
states in the formalism. Spectra of field free molecules, 
like that of copper trimer [3l or sodium trimer (l5| 
have been shown to bear the signature of the B erry 
phase in the electronic wave function due to CIs |16[ ■ 
The same classic consequence of the Berry phase in the 
electronic wave function is also found here for light in- 
duced CIs. In sharp contrast to field free polyatomic 
molecules where the CI is given by nature, the en- 
ergetic position of the CI discussed here can be con- 
trolled by the laser frequency and the strength of its 
non-adiabatic coupling by the laser intensity. Finally, 
we stress that for polyatomics there are, of course, also 
laser induced CIs. It is anticipated that their inter- 
play with the CIs given by nature will lead to a wealth 
of new phenomena. 
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